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We perform Monte Carlo study of the 3d jV — 6 superconformal U (N) x U (N) Chern-Simons 
gauge theory (ABJM theory), which is conjectured to be dual to M-theory or type IIA superstring 
theory on certain AdS backgrounds. Our approach is based on a localization method, which 
reduces the problem to the simulation of a simple matrix model. This enables us to circumvent the 
difficulties in the original theory such as the sign problem and the SUSY breaking on a lattice. The 
new approach opens up the possibility of probing the quantum aspects of M-theory and testing the 
AdS/[ I CFT3 duality at the quantum level. Here we calculate the free energy, and confirm the N 3 / 2 
scaling in the M-theory limit predicted from the gravity side. We also find that our results nicely 
interpolate the analytical formulae proposed previously in the M-theory and type IIA regimes.* 
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1. Introduction 

By now various regularization methods for supersymmetric gauge theories have been found, 
such as the lattice regularization (see e.g., ref. [jl|]), Fourier mode regularization [Q], the large-A^ 
reduction [|[ |J and non-commutative geometry However, all these methods require a lot of 
computational costs due to the existence of the dynamical fermions. In this article we introduce a 
new simulation method for investigating a class of supersymmetric field theories via localization 
method which reduces the evaluation of certain observables to calculations in simple matrix 
models. As a demonstration, we present numerical results [Q] for the so-called ABJM theory [j|], 
which is the 3d JV = 6 superconformal U (N) x U (N) Chern-Simons gauge theory. 

2. Localization method for general 3d jV = 2 supersymmetric field theory on S 3 

In this section we explain the basic idea of the localization method [||] and write down the par- 
tition function of a general 3d jV = 2 supersymmetric field theory on S 3 in terms of a matrix model 
[§]. The ABJM theory belongs to this class of theories. The localization method has been applied 
[|]] to 4d j¥ = 4 super Yang-Mills theory, and some conjecture on the half-BPS Wilson loops 1 
[11] has been confirmed. Those readers who are interested in just understanding our numerical 
results may skip this section. 

Let us consider the partition function of a supersymmetric field theory, 

f*r s[ * ] , (2.1) 

where <I> represents the collection of the components fields. Let us suppose that the action is 
invariant under an off-shell supercharge Q, namely 2 <2>S[ C I > ] = 0. Then, the closure of the SUSY 
algebra requires Q 2 = J^fg, where Jzfg is the generator of a bosonic symmetry the theory has. The 
first step of the localization method is to consider the deformation by a Q-exact term as 

Z(f) = J@<[>e- s W- tQV W, (2.2) 

where V is any fermionic functional satisfying J^VfO] = 0. By taking the derivative with respect 
to t, we obtain 



dZ(t) 
dt 



J ^ {QVme-^-WW = -J 9® Q (yme-W-WW} 

= J (Q®&) V [d^-sra-'G^*] . (2.3) 



If we assume the Q-invariance of the measure {Q3)<& = 0), namely that Q is non-anomalous, then 
the deformed partition function Z(t) should be independent of the parameter t. This implies that 
the original partition function Z can be written as 

Z = lim Z(t) = Z(t) = lim / e sm-tQvm . (2.4) 



1 This formula is also reproduced by a numerical simulation in the large-A' limit [|k 
2 Here we assume the absence of the boundary term. 
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In this limit, the saddle point approximation around the classical solution to QV = becomes exact. 
Hence we obtain 

Z = £e X p(-S[<J>o])Zi-ioop(<*>o), (2.5) 

where <J>o is the 'localized' configuration determined by ( 6 V )[*£>()] = 0. The summation £<£ over 
the saddle points should be understood as an integration if the saddle points are labeled by contin- 
uous parameters. The one-loop determinant Zi_i oop around <I>o is given by 



(2.6) 



Z,_ loop = lim f @{S<S>) e-WW 

We can also use this method to calculate ^-invariant operators such as supersymmetric Wilson 
loops |§| . 

Let us apply the localization method to a general 3d jV = 2 supersymmetric gauge theory on 
S 3 which is a Yang-Mills Chern-Simons gauge theory with arbitrary gauge group G = G\ x • • • x G r 
and Chern-Simons levels coupled to arbitrary number of .jV = 2 chiral multiplets with arbitrary 
representations and R-charge assignment 3 . The formula for the partition function is obtained as [pi] 



1 

W\ 



jrankG, a (l) jrankG,- (r) J_ 

- n A Vec(^ (fl) )n A Mat(^;^)> 



(2n 



,rankG 



(271 



,rankG r 



(2.7) 



0=1 



where \W\ is the order of the Weyl group of G, and (jW is the Cartan part of the adjoint scalar 
in the vectormultiplet with the gauge group G a at the localization point. Ay",(a( fl )) represents the 
contribution from the vector multiplet with the gauge group G a given by 4 



n 



2 sinh ■ 



a 



• exp 



_k<j(«) . CT («) 
4n 



(2.8) 



where labels the positive roots of G a . A^ t (a;q a ) is the contribution from the chiral multiplet 



with the representation 2% a and R-charge q a ( q a = 1/2 in the canonical assignment) 



2% 



(2.9) 



where p a is the weight vector of M a and f(z) is given by 



f(z) = exp 



-zdog(l 



„2izz 



(2.10) 



As a special case of a pair of chiral multiplets with the representation M and M in the canonical 



R-charge assignment, which coiTesponds to the jY = 4 hypermultiplet, the formula ( |2.9| ) reduces 
to the following simple form 

1 



Af at (a;l/2)Af at (a;l/2)= n 



? cosh 



(2.11) 



More generally, we can also include mass and FI terms [H], 
4 Note that this formula is independent of the Yang-Mills gauge coupling. This is because we can choose Q- V[<f>] as 
the action of JV = 2 super Yang-Mills theory itself. Then the deformation parameter t is nothing but the gauge coupling. 
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3. Numerical methods for the ABJM theory at arbitrary N and k 

Now let us consider the ABJM theory, which is the 3d jY = 6 superconformal U (N) x U (N) 
Chern-Simons gauge theory [gj. The Chern-Simons levels (the analogue of the gauge coupling 
constants) corresponding to two gauge groups are quantized to be integers, k and —k. This theory 
is conjectured to be dual to M-theory on AdSi, x S 1 jTL k for k <C N 1 ^ 5 , and to type IIA superstring 
on AdS^ x CP 3 at k -C N <C k 5 . The planar large-A^ limit is defined as the large-A^ limit with the 't 
Hooft coupling constant X =N/k kept fixed. 

The Monte Carlo study of the ABJM theory by usual lattice approach seems quite difficult for 
the following three reasons. Firstly, the construction of the Chern-Simons term on the lattice is 



not straightforward, although there is a proposal [12]. Secondly, the Chern-Simons term is purely 



imaginary in the Euclidean formulation, which causes the sign problem in the importance sampling. 
Thirdly, the lattice discretization necessarily breaks supersymmetry, and one needs to restore it in 
the continuum limit by fine-tuning parameters 5 . 

In order to circumvent these problems, we apply the Monte Carlo method to a matrix model 



obtained via the localization. According to the general formula ([277 ), the partition function of the 

Vi-Vj 



ABJM theory on S 3 is given by 2 



1 f d N n d N v rL<y 

Z(N,k) - 



2 sinh 



2sinh ■ 



(N\) 2 J (2ti) n (2k) n 

1 lij 



e&HLrtf-yn, (3 .i) 



2 cosh 1 



which is commonly referred to as the ABJM matrix model. From the partition function, we define 
the free energy as 

F{N,k) =\ogZ(N,k). (3.2) 

Thus the ABJM free energy is given just by a 27Y-dimensional integral. Note that the ABJM matrix 
model describes the continuum theory without any regularization artifact. 



The ABJM matrix model in the form (3.1) is not suitable for Monte Carlo simulation since the 



integrand is not real positive. However, as we reviewed in Appendix B of [Q] in detail (See also the 
original work [|14|]), one can rewrite the ABJM matrix model as follows. 

Z(N,k) = C N , kg (N,k), C^ = — L^, 



FLwtanh 2 ( %^ 

>(N,k) = I d N x *\ u ) U J = I d N xe s(N^ u -^)_ (3 .3) 



fT,-2cosh(x i -/2) 

An important point here is that, in this form ( |3.3[ ), the integrand is real positive, and we can perform 
Monte Carlo simulation in a straightforward manner. 

In order to calculate the partition function, we need to rewrite it in terms of expectation values 
of some quantities, which are directly calculable by Markov-chain Monte Carlo methods. The 
basic idea is to calculate the ratios of the partition functions for different N as expectation values 6 . 
Let us decompose ,/V into N = N\ + N2 and consider the ratio 



5 This might be overcome by a non-lattice regularization of the ABJM theory Jl3| ] based on the large-iV reduction 
on S 3 which is shown to be useful in studying the planar limit of the 4d .jV = 4 super Yang-Mills theory |Kjj. 
6 We can also calculate the ratios of the partition functions for different k as expectation values ]m. 
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Figure 1: (Left) The free energy is plotted against N 3 I 2 for k = 1,2,4,6,8. The data points can be fitted 
to straight lines, which implies F ~ N 3 / 2 as N increases. (Right) The M-theory limit of the free energy 
limm^ooF /N 3 / 2 is plotted against \f~k~. Our data are in good agreement with the result (5.1) predicted from 
the eleven-dimensional supergravity, which is represented by the solid line. 



i(N,k) 



^(Nuk-xi,- ,x Nl )+S(N 2 ,t,x Nl+ i,- ,x N )-S(N,k) \ (2.4) 



g(N u k)g(N 2 ,k) /n u N 2 

x 7 / Ni ,N 2 

where the symbol (• ■ ■ )tf l jf 2 represents the expectation value with respect to the "action" given by 




S(Ni,k;xi, ■ ■ ■ ,Xnj ) + S(N2,k;xN l+ i, ■ ■■ ,xn). In order to calculate the right-hand side of (3.5) with 
good accuracy, it is necessary to take N 2 small enough to make sure that the observable in ( |3"3| ) 
does not fluctuate violently during the simulation. In actual calculation we use N 2 = 1. Then, by 



calculating (3.5) for Ni = 1,2,3, •• • and by using the N = I result g(l,k) = K, we can obtain the 



free energy for iV = 2,3,4, • • ■ successively with a fixed value of k. 

4. Results for the free energy 

We present our numerical result [Q] for the free energy of the ABJM theory. First we consider 
the large -N limit with fixed k, which is conjectured to correspond to the eleven dimensional super- 
gravity on A0IS4 x 5 7 /Zjt. In refs. [|15|, [16|, 17], the free energy in the M-theory limit (N — > 00 with 



k fixed) has been calculated by various analytic methods and confirmed the prediction 

x\/2k 



F SU GRA = ^— W 3/2 (4.1) 

from the dual eleven-dimensional supergravity. Figure [I] (Left) shows that the free energy F grows 
in magnitude as 7V 3/2 with N. Actually F /N 3/1 behaves as F (N, k) /N 3/2 = h (k)+h x (k) /N, which 
enables us to obtain the M-theory limit ho(k) = lim^oo/ 7 '(N ,k) / 'N 3 '/ 2 reliably. In fig. [j] (Right) we 
plot ho (k) against yk, which confirms the prediction ( |4.l| ) from the eleven-dimensional supergrav- 
ity for k = 1 , 2, • • • ,10 very precisely. 

Let us next study the finite-TV effects. An important analytical result on finite N effects is that 



the 1 /N corrections around the planar limit are resummed in a closed form [Nq, 16] 
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Figure 2: (Left) The free energy of the ABJM theory for N = 4 is plotted against 1 /vA. The solid line 
represents the FHM result. The dotted line represent the perturbative results /\y e ak = 

2 log G% (N + 1) with the Barnes G-function G2 (x) . (Right) The difference F — Ffhm is plotted against N for 
various values of k. It reveals non-negligible discrepancies for each k, which are almost independent of N. 



Ffum(N,X) = log 



1 (4n 2 N 
71 



1/3 



Ai 



kN 2 

7^ 



2/3 



1 

24 



3N 2 



(4.2) 



where Ai(x) is the Airy function and the type of correction 0(e 2n ^-) is neglected. In fig. || 
(Left) we plot our results for = 4 and compare them with the FHM result (JO). We find that 



our result agrees reasonably well with the FHM result in the strong coupling regime. To see more 
precisely, we plot in fig. |2| (Right) the difference between our result and the FHM result against 
Af for various k. It turns out that there are discrepancies which are almost independent of N. 
This strongly suggests that the FHM result correctly incorporates the finite N effects except for a 
term which depends only on k. Note that this discrepancy cannot be explained by the worldsheet 
instanton effect 0(e~ 27C ^), which is neglected in FHM. See ref. [Q] for a natural interpretation of 
this discrepancy from topological string theory. 



5. Summary and discussions 

In this paper we have established a simple numerical method for studying the ABJM theory 
on a three sphere for arbitrary rank N and arbitrary Chern-Simons level k. The crucial point is that 
we are able to rewrite the ABJM matrix model, which is obtained after applying the localization 
technique, in such a way that the integrand becomes positive definite. By using this method, we 
have confirmed from first principles that the free energy in the M-theory limit grows proportionally 
to Af 3 / 2 as predicted from the eleven-dimensional supergravity. We have also found that the FHM 
formula with the additional terms describes the free energy of the ABJM theory in the type II A 
superstring and M-theory regimes. While we have focused on the free energy as the most funda- 
mental quantity in the ABJM theory, our method can be used to calculate the expectation values of 
BPS operators. For instance, it is possible to calculate the expectation value of the circular Wilson 



loop for various representations [19]. 
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We hope that the results of this work are convincing enough to show the power of the combi- 
nation of the localization method and numerical simulation. We expect further numerical study of 
various localized matrix models will reveal exciting new aspects of supersymmetric gauge theories 
and quantum gravity. 
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